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We investigate hole spin relaxation in intrinsic and p-type bulk GaAs from a fully microscopic 
kinetic spin Bloch equation approach. In contrast to the previous study on hole spin dynamics, we 
explicitly include the intraband coherence and the nonpolar hole-optical-phonon interaction, both 
of which are demonstrated to be of great importance to the hole spin relaxation. The relative 
contributions of the D'yakonov-Perel' and Elliott- Yafet mechanisms on hole spin relaxation are also 
analyzed. In our calculation, the screening constant, playing an important role in the hole spin 
relaxation, is treated with the random phase approximation. In intrinsic GaAs, our result shows 
good agreement with the experiment data at room temperature, where the hole spin relaxation is 
demonstrated to be dominated by the Elliott- Yafet mechanism. We also find that the hole spin 
relaxation strongly depends on the temperature and predict a valley in the density dependence of 
the hole spin relaxation time at low temperature due to the hole-electron scattering. In p-type GaAs, 
we predict a peak in the spin relaxation time against the hole density at low temperature, which 
originates from the distinct behaviors of the screening in the degenerate and nondegenerate regimes. 
The competition between the screening and the momentum exchange during scattering events can 
also lead to a valley in the density dependence of the hole spin relaxation time in the low density 
regime. At high temperature, the effect of the screening is suppressed due to the small screening 
constant. Moreover, we predict a nonmonotonic dependence of the hole spin relaxation time on 
temperature associated with the screening together with the hole-phonon scattering. Finally, we 
find that the D'yakonov-Perel' mechanism can markedly contribute to the hole spin relaxation in 
the low density case at moderate temperature and in the high density case at low temperature, 
where the Elliott- Yafet mechanism is suppressed due to the relatively weak scattering. 



PACS numbers: 72.25.Rb, 71.70.Ej, 71.10.-w, 71.55.Eq 



I. INTRODUCTION 

In the past decade, semiconductor spintronics, with 
the aim of reahzing favourable devices for future applica- 
tion based on the spin degree of freedom, has attracted 
much attentionjir— One of the most critical challenges for 
such devices hes in the control of the spin lifetime, which 
is hmited by the unavoidable spin relaxation and/or de- 
phasing process in semiconductors. In this sense, under- 
standing of the carrier spin relaxations and/or dephasings 
is a critical issuef2^" In bulk materials, the spin relax- 
ation and/or dephasing properties of the electrons and 
the underlying physics have been well understood after a 
long-time researchi^^— However, the study on hole spin 
dynamics, which occurs on the time scale of the momen- 
tum scattering time (usually < 1 ps), is fairly rare, partly 
because of the limited resolution on the detection of such 
an ultrafast process. To our best knowledge, the hole 
spin lifetime in bulk semiconductors was only measured 
in intrinsic GaAs at room temperature, which evaluated 
the spin relaxation time of the heavy hole (~ 110 fs).^^ 

One picture to explain the picosecond hole spin relax- 
ation (HSR) time is the D'yakonov-Perel'-like description 
associated with the momentum scattering and the spin 
precession between the light-hole (LH) and heavy-hole 
(HH).^^ However, this picture works only if the broad- 
ening of the energy spectrum due to scattering is larger 
than the interhand splitting, corresponding to the vicin- 



ity of the zone center. Otherwise, holes are driven to 
large momentum states where the degeneracy between 
the LH and HH bands is significantly lifted, hence it is 
more appropriate to treat the LH and HH bands sepa- 
rately. The HSR of the individual LH and HH bands 
was studied by Yu et al.^ who obtained the accurate 
band structure from the tight-binding Hamiltonian with 
the spin-orbit coupling (SOC) and calculated the spin 
relaxation time of the LH (HH) as the inversion of the 
decay rate of the quasi-spin polarization, that is, the 
population difference of the LH (HH) between the two 
quasi-spin bands. The Elliott- Yafet (EY) mechanism^ 
there described as the direct quasi-spin-flip scattering, 
was claimed to be the solo mechanism for the HH and 
the dominant one for the LH, while the D'yakonov-Perel' 
(DP) mechanism^ treated as the intraband precession 
together with the quasi-spin-conserving scattering from 
the motional narrowing relation was found to be unim- 
portant. Since the calculation was based on the single 
particle theory, the Coulomb interaction which has been 
demonstrated to result in intriguing many-body effects 
during electron spin dynamics^i^^— was missed in that 
work. Moreover, the HSR time was extracted from the 
quasi-spin polarization instead of the exact spin signal. 
The feasibility of this approach needs to be verified. In 
other words, the results obtained there is not exactly 
the HSR time, but the quasi-spin relaxation time. Re- 
cently, the hole spin dynamics in intrinsic bulk GaAs 



with the Coulomb interaction was microscopically inves- 
tigated from an eight-band Kane Hamiltonian by Kraufi 
et al-r^ where the HSR time, directly from the decay of 
the hole spin expectation, was shown to be quite close to 
the experimental result. It was also shown that the HSR 
time can be slightly different from the quasi-spin relax- 
ation time. Thus, it seems that ultrafast HSR in bulk 
zinc-blende semiconductors has been successfully inter- 
preted theoretically. However, one may notice that the 
intraband coherence, i.e., the non-diagonal components 
of the density matrices was missing in that work, which 
actually will lead to two consequences: the exaggeration 
of the EY mechanism and the exclusion of the DP mech- 
anism. Moreover, the nonpolar interaction of holes with 
transverse- and longitudinal-optical-phonons due to the 
deformation potential coupling, which was shown to pro- 
duce a significant contribution on the charge dynamics 
of holes, •^^"^^ was not included. In this sense, the results 
in Ref. 26 should be reexamined. 

In the present work, we employ the fully microscopic 
kinetic spin Bloch equations (KSBEs)^'i^>2^'^>^ to in- 
vestigate the HSR due to the EY and DP mechanisms 
in intrinsic and p-type bulk GaAs. The EY mechanism 
here corresponds to the decay of the spin polarization 
solely due to the scattering associated with the inter- 
band mixing, between not only the HH and LH bands 
but also the conduction and valence bands; whereas the 
DP mechanism stands for the additional contribution by 
including the intraband spin precessions. We first ana- 
lytically derive the valence band structure from the four- 
band Luttinger HamiltonianH together with the Dressel- 
haus SOC (Ref.[33) due to the bulk inversion asymmetry 
in zinc-blende semiconductors. We demonstrate that the 
Dresselhaus SOC induces an intraband splitting between 
the two HH bands as well as the two LH bands with cu- 
bic wave- vector dependence, which supplies an additional 
HSR channel due to the DP mechanism. We find that 
the intraband splitting of the HH bands vanishes under 
the spherical approximation of the Luttinger Hamilto- 
nian, and then the DP mechanism becomes irrelevant 
to the HH spin relaxation.^'' This implies the limitation 
of such an approximation in calculating the HSR time. 
Therefore, we later obtain the intraband splitting and 
wave functions by diagonalizing the full eight-band k • p 
Hamiltonian beyond the spherical approximation. Then 
we calculate the HSR time from the KSBEs with the in- 
traband coherence and all the relevant scatterings, such as 
the hole-impurity, hole-hole, hole-electron, hole-acoustic- 
phonon, polar and nonpolar hole-optical-phonon scatter- 
ings, explicitly included. In the intrinsic case at room 
temperature, the HSR is dominated by the EY mech- 
anism and the HSR time shows good agreement with 
the experiment. We find that the HSR time can be sig- 
nificantly manipulated by changing the temperature and 
excitation density. In the p-type materials, we predict in- 
triguing nonmonotonic behaviors of the HSR time in both 
the density and temperature dependences. We show that 
the nonmonotonic features reflect the role of the screen- 



ing. Moreover, we find that the EY mechanism is usually 
the major mechanism of the HSR, but the DP mecha- 
nism can still be comparable with the EY mechanism in 
some special cases, such as the high density regime at 
low temperature and the low density regime at moderate 
temperature. 

This paper is organized as follows. In Sec. II, we set up 
our model and derive the effective Dresselhaus field from 
the Luttinger Hamiltonian. The KSBEs are also con- 
structed in this section. In Sec. HI, we investigate the 
HSR in both intrinsic and p-type GaAs. The compari- 
son of the calculations with and without the intraband 
coherence is given in this section to illustrate the role of 
the intraband coherence. The relative contributions of 
the DP and EY mechanisms are also discussed. Finally, 
we summarize in Sec. IV. 



II. MODEL AND KSBE 

To qualitatively analyze the band structure and the 
intraband splitting of the HH and LH bands, we start 
from the perturbation method with the 4x4 Hamiltonian 
(in the basis of the eigenstates of Jz with eigenvalues 
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The first term on the right-hand side of the equation 
is the Luttinger Hamiltonian,"^^ where F — —^^^[{li + 
I2){kl + kl) + (71 - 2^^)k% G = ^£-[{^,-Z){kl + 
fc2) + (71 + 272)fc2], H = 2V32^73(fcx - iky)k,, and 

^ = 2^[^72(fc? - kl) - i2^/i-^■ik^ky\. 7, are Kohn- 
Luttinger parameters. We denote this term as H^ in the 
following. The second term is the Dresselhaus SOC of the 
valence band, H^y^^ = hk • J, with the effective magnetic 
field^^^ 



hk — O41 [kx{ky — k^),ky[k^ — k^),kz(k^ ~ ky)\. (2) 

Ji represent the spin-3/2 angular- momentum matrices. 
It is obvious that the band structure is mainly determined 
by Hq in the vicinity of the zone-center. The energy 
spectrum from Hq readsSIiiS, 

Ek/iM - -2&{7ifc'T[47|fc* + 12(7|-7l) 

x{klkl + klkl + klkl)]'2}, (3) 

and the wave functions can be expressed by 
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By identifying ip-^ and '02 as pseudo-spin- up and 
-down states separately, one obtains the effective mag- 
netic field of the Dresselhaus SOC {^h/i) between them. 
Specifically, one rewrites HgvSv in the representation of 
{ipi }, and derives individual HH and LH blocks from 
the Lowdin partitioning method^ upto the cubic power 
of the wave vector, 

Hh/iik) = Eh/iMhx2 + ni^iik) ■ (T, (6) 

where <t are the Pauli matrices. The expressions of 
the effective magnetic field together with the coefficients 



h/l 



6j and c^ are given in AppendixlX] The pres- 



^h/lM 



ence of the intraband splitting AEfi/i,k 
^h/i]i ~ 2|r2^/;(k)| obviously indicates that the DP 
mechanism can contribute to the HSR in the presence 
of the scatterings 1^*22, Interestingly, one finds Q.\ van- 
ishes once the spherical approximation (corresponding to 
72 = 73) 38,39 jg applied, which means that the contribu- 
tion of the DP mechanism to the HH spin relaxation is 
lost within the spherical approximation scheme. In other 
words, the difference between 72 and 73, arising from the 
remote bandsf^i is important in counting the HH spin 
relaxation. 

We should point out that the above perturbation ap- 
proach gives the precise value of the intraband splitting 
only in the vicinity of the zone-center. For the regime far 
away from the center, the modification from the split- 
off and conduction bands should be considered. There- 
fore, the intraband splittings Ai<J^ « ^ are obtained from 
the diagonalization of the 8x8 Kane Hamiltonian Hk 
(Ref..40.) in our numerical calculation, that is, by solving 
the Schrodinger equation 
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with i — 1,2. 
sity matrices 



representation,— i.e 



/j/;,k/ — ^/i//,klS/i//,k/' 

We define the single particle den- 
as 4 x 4 matrices in the helix 
under the basis of the eigenstates 



llCh.k)' l?^,k)' IC/,k)' ICi\)}; The unitary transformation 
from the helix representation to the collinear representa- 
tion, the basis of which is defined as the eigenstates of 
the angular momentum operator J^, is given by 

Pk = V^pIvI (8) 

with C/u = (eU,e^,k,e/,k,^F,k)- The KSBEs in the helix 
representation from the nonequilibrium Green-function 
method reads^i^^ 



dtpl = dtpl\^^d,pl\ 
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The coherent term can be written as 
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In the helix frame, ff''(k) = diag(ii;i_k,i5^,k,^',k>£^f,k): 
and Eq. (fTO)l then can be rewritten as dtp\(m^n)\ , — 
~ip^{m,n){Em k — En k)- The scattering term is given 
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l~Pk- Tll{m,n) = Sm,r^Sn,' 



where, p^^ = p'^ and p^^ 

The hole- impurity scattering matrix element jt/q, — 
Z^V^ with Z taken to be 1 in the calculation. V^ = 
e^/[Ko((7^ + K^)]. kq denotes the static dielectric con- 
stant. Here, the screening constant k is calculated from 
the random phase approximation (RPA)ji^ The detail 
of the polar carrier-longitudinal-optical(LO)-phonon and 
carrier-acoustic(AC)-phonon scattering elements |Mq_Ap 
can be found in Refs.124! and \^. Besides, the longitu- 
dinal and transverse optical modes can also contribute 
to the nonpolar hole-optical-phonon scattering. The ma- 
trix elements of these scatterings are givenSI by M^°"'' — 

( 



"q,A 



-yf^Dq.x, where D in the square root represents 



the crystal density. The potential matrix Dq\ is given 



in Appendix[Bl N^^ 



[exp{ujq^x/kBT) - ly 
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The function Q in the Coulomb scattering term reads 
Qiq,w) = ^ (5(^,,k" - Er,'k"-ci - w) 

777;' k" 
XTr[S'k"-q,k"Pk"^T^"5'k",k"-qrk'„_qPk;,1q] 

+ Y. ^i^-^" - ^-k"-q - u;)Tr[p>k„p<k._q].(12) 

k" 

The second term on the right-hand side of the above 
equation describes the contribution of the hole-electron 
scattering with pe.k representing the electron density ma- 
trices. 



III. NUMERICAL RESULTS 



with [, ] representing the commutator. Sk'.k = Ul^iUk 
The first term on the right-hand side of the equa- 
tion comes from the Coulomb Hartree-Fock contribu- 
tion, which can be neglected for small spin polarizatioui^ computation are all taken from Ref.l44l. By numeri 



In this section, we present our results in bulk GaAs 
with the measured value of the optical deformation po- 
tential do = 48 eVi2Si^ The other parametersii in our 
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FIG. 1: (Color online) (a) Energy spectra and (b) intraband 
splittings of the valence bands (solid curves) along the [110]- 
direction from the diagonalization of the 8x8 Kane Hamil- 
tonian. The dashed curves in (b) show the results from the 
4x4 Luttinger Hamiltonian. a is the lattice constant, (c) and 
(d) show the intraband splittings of the HH and LH bands at 
|k| = 0.l7r/a as function of the wave direction {0,tp). 



the effective-mass approximation. Moreover, the split- 
ting between the HH and LH bands are tens of meV even 
for the Fermi energy smaller than 10 meV, corresponding 
to n/i « 3 X 10"'^^ cm~^ at K. Therefore, we neglect the 
interband coherence between the HH and LH bands and 
reduce the 4x4 hole density matrices into the LH and 
HH blocks, both are 2x2 matrices, in our computation. — 
In Fig.[TI^b), we illustrate the intraband splittings of the 
three valence bands, i.e., A_E^^k = -^-^ k ~ -^77 k' from the 
8x8 Kane Hamiltonian along the [110]-direction as solid 
curves. Here, E'i (E^) is the larger (smaller) eigenvalue 
of the r/-band with 77 corresponding to the HH, LH and 
split-off bands. One can see that the splitting of the HH 
band is much smaller than that of the LH band. This can 
be understood as follows: Within the spherical approx- 
imation, the LH band itself presents an intraband split- 
ting due to the Dresselhaus SOC, while the HH band is 
still doubly degenerate as mentioned in Sec. II. However, 
the anisotropy property of the valence band makes the 
HH states contain some LH components so as to lift the 
degeneracy of the HH band. In other words, the SOC 
induces the splitting of the HH band indirectly, hence 
the magnitude is smaller than the direct splitting of the 
LH band. For comparison, we also plot the intraband 
splittings from the perturbation approach up to the or- 
der of k^ based on the 4x4 Luttinger Hamiltonian in 
Fig.[TJb) as dashed curves. One can see that the pertur- 
bation approach only performs well for k < 0.167r/a. For 
larger wave vectors, the cubic wave- vector dependence of 
the intraband splitting is violated due to the interband 
coupling. The intraband splittings of the HH and LH 
bands are plotted as function of the direction of the wave 
vector in Fig.[ljc) and (d), respectively, from which the 
anisotropy properties can be clearly seen. 

Since the intraband splitting is much smaller in the 
density and temperature regimes studied in the present 
work compared to the Fermi energy, one can neglect 
the intraband splitting in the 5- function in Eqs. (fTTj) 
and p^ ."^^ For the numerical treatment of the scatter- 
ing term, we take the isotropic energy spectrum from 
the effective-mass approximation, i?^ y^ j, « ^h/i k ~ 

Eh/lm = -2&(7i T 7)fc' with 7 = (72 + 73)/2^^ The 
feasibility of this widely adopted approximation has been 
shown in the literature. ^^"'*^ Moreover, we find that the 
density of states from this approximation is almost the 
same as that obtained from the anisotropic spectrum. 
The average momentum scattering time in this scheme is 
also comparable to that from the real band structure. 



cally diagonalizing the 8x8 Kane Hamiltonian, one ob- 
tains the energy spectra of the valence bands. The re- 
sults along the [110]-direction are plotted in Fig.[T][a), 
where the HH, LH and split-off bands all present intra- 
band splittings. One notices that both the HH and LH 
bands are approximately parabolic within 0.15 eV from 
the valence-band top, which indicates the feasibility of 



A. HSR in intrinsic GaAs 

In this part, we investigate the HSR in intrinsic case. 
Our discussion is based on two physical quantities, i.e., 
the quasi-spin polarization and the spin polarization. 
The former describes the population difference between 



the two HH and LH quasi-spin bands, defined as 

^^ = E(ftV(l'l)-/^'k(2,2))/n,, (13) 

77k 

with nil being the hole density. The latter is calculated 
as the spin polarization along the [001]-direction 






(14) 



which reflects the optical orientation signal in experi- 
ment. 
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FIG. 2: (Color online) Temporal evolution of the normal- 
ized quasi-spin polarization (dashed curves) and spin polar- 
ization (solid ones) at room temperature with the hole den- 
sity Uh = 10^^ cm"'' in intrinsic GaAs. The curves labelled 
as "full" ("incoherent") are obtained from the full calcula- 
tion with (without) the intraband coherence, while the ones 
with "EY" show the results without the intraband preces- 
sion term, corresponding to the contribution from the EY 
mechanism. The results from the calculation including the 
intraband coherence but without the nonpolar hole-optical- 
phonon scattering are plotted as the curves denoted as "with- 
out NPHOP" . The dotted curve represents the spin polariza- 
tion from the calculation scheme used in Ref. 26, i.e., without 
the nonpolar hole-optical-phonon scattering and intraband 
coherence. The relaxation times are given correspondingly. 

In our computation, we take the initial state from the 
optical orientation due to the pump pulse with the small 
polarization P°p* = (/+ - /")/(/+ + /") = 2 %,iS where 
/"*" and /~ are the intensities of the a^- and cr"-polarized 
light. In that case, one sets the electron density matri- 
ces as the Fermi distribution with the spin polarization 
1 % at the lattice temperature and keeps it unchanged, by 
taking into account of the fact that the electron spin (mo- 
mentum) relaxation time is much longer (shorter) than 
the time scale of the HSR. The initial hole density ma- 
trices are also set to obey the Fermi distribution in the 
collinear spin space to describe the optically excited con- 
dition. By solving the KSBEs, one obtains the tempo- 
ral evolution of the spin polarization P* (quasi-spin one 
P^), from which the spin (quasi-spin) relaxation time is 
extracted. 



We first take the hole density n^ = 10""^^ cm~^. Fig- 
ureO illustrates the temporal evolution of the spin po- 
larization (normalized by the value at i = 0) from the 
full calculation at room temperature as the red solid 
curve (labelled as "full"). One can see that the spin 
polarization decays exponentially with the HSR time 
Ttot ~ 0.11 ps, which agrees perfectly well with the 
experimental value of the HSR time, O.llilO % ps.— 
By removing the coherent term from the KSBEs, one 
switches off the DP mechanism and obtains the results 
solely due to the EY mechanism as the blue solid curve. 
The HSR time tey in this case is about 0.13 ps. Since 
the EY mechanism is irrelevant to the spin precession. 



from 



'EY 
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one extracts the spin lifetime 



due to the DP mechanism tdp ~ 0.72 ps, which is much 
longer than that of the EY mechanism. Therefore, the 
EY mechanism is dominant in this case. Moreover, we 
find that the relaxation time of the quasi-spin polariza- 
tion (red dashed curve labelled with "full" ) is very close 
to that of the spin polarization. 
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FIG. 3: (Color online) HSR time from the full calculation as 
function of the hole density Uh in intrinsic GaAs. The result 
from the calculation without the hole-electron scattering at 
50 K is labeled as "without h-e". The dot represents the ex- 
periment data at room temperature with the photo-excitation 
density n^ = 10^^ cm~^ fRefflal. In the calculation, the elec- 
tron density Tie = Uh. 

In the figure, we also plot the results from the calcu- 
lation without the intraband coherence as orange curves 
(labelled as "incoherent"), where the non-diagonal ele- 
ments p,,k(l,2) are artificially set to be zero during the 
calculation as done in Ref.[2a. Compared with the results 
labelled with "EY" , one immediately finds that this pro- 
cedure markedly exaggerates the contribution of the EY 
mechanism, (The DP mechanism vanishes in this scheme 
obviously.) which clearly shows the problem of the miss- 
ing of the intraband coherence. The reason lies in the fact 
that the spin-conserving scattering becomes the spin-flip 
one once the non-diagonal elements of the density ma- 
trices in helix representation are neglected4i In Ref.[26l. 
the authors argued the disregarding of the intraband co- 
herence by arguing that the optically driven coherence 



between single particle states would be very small and 
the SOC could not drive any coherence further during 
the time evolution of the density matrices in the helix 
representation (i.e., under the basis called "intelligent 
basis" there). However, this is actually incorrect. Since 
the SOC information is transferred into the wave func- 
tions of the eigenstates in the helix representation, any 
spin-conserving scattering events can induce the intra- 
band and/or interband coherence4i Therefore, the same 
HSR time is obtained if one chooses the helix initial con- 
dition, corresponding to the initialization without any 
non-diagonal elements of the density matrices a.t t — 0. 
In the remaining part of this paper, all the results are 
obtained from the calculation including the intraband co- 
herence. 

The green curves (labelled as "without NPHOP") in 
Fig.[2]obtained from the calculation without the nonpolar 
hole-optical-phonon scattering give the HSR time around 
0.17 ps, much longer than that from the full calculation. 
This clearly shows the importance of the nonpolar hole- 
optical-phonon scattering in the HSR. It is noted that 
this scattering is also missed in Ref.l26l For comparison, 
we further plot the spin polarization from the calculation 
without the nonpolar hole-optical-phonon scattering and 
the intraband coherence as black dotted curve, which cor- 
responds to the calculation scheme in Ref.[2a. It is seen 
that the missed terms have marked contribution to the 
HSR. 

Since the direct relation between the EY mechanism 
and the scattering strength, the HSR time due to the 
EY mechanism decreases (increases) when the scatter- 
ing becomes stronger (weaker). To show the influence 
of the scattering more clearly, we vary the density and 
temperature. The density dependence of the HSR time 
is plotted in Fig. [3] One can see that the HSR time at 
room temperature monotonically decreases with increas- 
ing photo-excitation density. Since the electrons and 
holes are both in the nondegenerate regime, the Coulomb 
scattering rate becomes larger as the density increases, 
according to the estimation of the average Coulomb scat- 
tering rate (r^^)"^ ex n/^/T^/^^ As a result, the HSR 
time decreases. One can also see that the HSR time 
markedly increases with decreasing temperature, due to 
the suppression of the hole-phonon scattering. 

More interestingly, we predict a valley in the density 
dependence of the HSR time at low temperature. The 



minimum occurs at Uh 



= 4 X 10^ 



cm 



(cor- 



responding to the Fermi temperatures Tp = 74 K for 
electrons and T^ = 9 K for holes) at T = 50 K and 
nh = ne = 2x IQi^ cm~3 (T| = 216 K and T^ = 27 K) 
at T = 100 K. The agreement between Tp at the val- 
ley and the lattice temperature T reveals that the valley 
results from the different density dependence of the hole- 
electron scattering in the degenerate and nondegenerate 
limits. As mentioned above, the HSR time decreases with 
increasing density in the nondegenerate limit. However, 
for high densities, electrons first enter into the degenerate 
regime. Therefore the hole-electron scattering strength is 



significantly decreased due to the Pauli-blocking of elec- 
trons. As a result, the HSR time due to the EY mech- 
anism increases in this regime, and finally results in the 
valley at the crossover between the degenerate and non- 
degenerate regimes of electrons. In the figure, we also 
plot the result from the calculation without the hole- 
electron scattering at 50 K. It is clear to see that the 
valley at n^ = rie = 4 x 10^^ cm~'^ disappears, as ex- 
pected. Instead, one finds that another valley shows up 
a,t Uh = Ue = 2 X 10^'^ cm~3, where Tp (27 K) is com- 
parable to the lattice temperature T (50 K). This means 
that holes also enter the degenerate regime at such a high 
density and the hole-hole scattering rate is then given by 



2/3 50 



Therefore, the HSR time increases 



(t,^)-i (X rVn 
with increasing the density for Uh = n^ > 2 x 10^^ cm^'^. 
One notices that the Coulomb scattering can also induce 
a peak (instead of valley) in the density dependence of 
the electron spin relaxation time limited by the DP mech- 
anism, first predicted by Jiang and Wui^ and realized 
experimentally by Kraufi et al.J^i^ 



B. HSR in p-type GaAs 

In this part, we study the HSR in p-type GaAs. We 
still apply the collinear initial spin polarization P* — 
2 %. As the photo-excited carrier density is much smaller 
compared to the doping density, we take Uh = Ui and 
neglect the hole-electron scattering in the computation. 



1. Density dependence 

We start from the density dependence of the HSR 
time. The results from the full calculation are plotted 
in Fig. Ufa), which shows rich and intriguing nonmono- 
tonic features. At low temperature, e.g., at 5 K, the 
HSR time first increases and then decreases with in- 
creasing hole density. As a result, a peak appears at 
Ufi — 4: X 10^® cm"-^, where the corresponding Fermi 
temperature T^ = 9 K is comparable to the lattice tem- 
perature. As the temperature increases, the peak moves 
to the higher density regime (n/j = 3 x 10^^ cm~^ with 
T^ = 35 K for 20 K and rih = W^^ cm^^ with T^ = 79 K 
for 50 K). Simultaneously, a valley gradually appears in 
the low density regime. Although the peak also locates 
at the crossover between the degenerate and nondegen- 
erate regime, the underlying physics of the peak is again 
different from the one in the electron spin relaxation time 
in intrinsic and n-type materials. -'^^ To explain the pres- 
ence of the peak, we turn to analyze the change of the 
scattering strength. Since the intraband splitting is too 
small to affect the HSR for the hole states near the zone- 
center, the HSR process is mainly determined by the EY 
mechanism instead of the DP mechanism for low density 
case, especially at low temperature. 

As a qualitative analysis for the case rii — n^, one 
can focus on the major scattering mechanism, i.e., the 




r^/^. For K ':§> q, one neglects the q^ term in r'" and 
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FIG. 4: (Color online) (a) HSR time from the full calculation 
as function of the hole density Uh in p-type GaAs. We take 
the impurity density rii = nn- The hole- impurity scattering 
time Tp ' at 20 K is plotted as the black dotted curve without 
symbol, (b) The ratio of the screening constant to the Fermi 
wave vector, K/kp, as function of the hole density. The solid 
line illustrates the level k = \/2fcF. 



hole-impurity scattering with the strength estimated by 
(t«)-i ~ m,^nh{kV^)H3ir^h^) ex nh{k{q^ + k^)-^), 
where k comes from the density of states and (...) stands 
for the average over the distribution. TOcs represents the 
effective mass and q is the momentum exchange. In the 
degenerate regime, only holes on the Fermi surface con- 
tribute to the HSR, therefore, one can estimate k ^ kp 
and q^ « 2k p with kp representing the Fermi wave vec- 
tor. Thus, one has 



r^' (X [2 



(K/kf 



(15) 



FigureHl^b) shows the ratio of the RPA screening con- 
stant K to kp. One finds that n/kp decreases with in- 
creasing density in the high density regime and r^ ~ t^ * 
decreases also. The density dependence of n/kp can 
be easily understood once the Thomas-Fermi screening 
(k cx n,/ ) is applied in the degenerate limit, ^^ which 

leads to n/kp cx n^ . 

However, the situation is quite different in the nonde- 
generate regime. In this limit, one has k cx n^ T^i/2 
according to the Debye-Hiickel screeningj^^ and q ^ k ^ 



obtains 



oc {nh{kn-^))-^ oc nnT-'''^, 



(16) 



which indicates that the HSR time increases with increas- 
ing density in this case. For k <^ q, the screening con- 
stant is neglected, and Tp* can be written as 



Tp^ CX {nh {kq 



'))- 



'/i 



l2i3/2 



(17) 



which decreases with increasing density. In Fig.|31Ja), we 
also plot the density dependence of Tp" at 20 K (the dot- 
ted curve without symbol), which agrees well with our 
discussion. It is clear to see that r^' is in the same order 
of magnitude as the HSR time t^ as expected. 

Now, the peak and valley in the HSR time in Fig.|3l^a) 
can be well understood. For example at 20 K, holes lie 
in the nondegenerate regime and the screening constant 
is small in the low density regime. As the density in- 
creases, the HSR time decreases according to Eq. (fT7)l . 
Nevertheless, with further increase of the density, k can 
be larger than q. Then the HSR time increases with den- 
sity according to Eq. p^ . By further increasing density, 
the system enters into the degenerate regime, and the 
HSR time decreases again following Eq. p3|) . By com- 
paring Fig. 3] (a) and (b), we find the crossover between 
Eqs. (HH]) and (|17p can still be qualitatively estimated by 
taking q^ ^ 2k p. At low temperature (5 K), the screen- 
ing constant is always large in the nondegenerate regime 
of our investigation [see Fig.|4jb)], so the valley is invis- 
ible in Fig.l^a). However, at high temperature (200 K), 
the screening is weak in both the degenerate and non- 
degenerate regimes and the HSR time monotonically de- 
creases with increasing density. 

To elucidate the role of the screening more clearly, 
we plot the results from different screenings at 20 K in 
Fig.EJa) . The green dotted and blue dashed curves corre- 
spond to the results with the Thomas-Fermi and Debye- 
Hiickel screenings, respectively. The total HSR time from 
the RPA screening is also plotted as the red solid curve. 
One can see that the result from the RPA screening 
is consistent with that from the Thomas-Fermi (Debye- 
Hiickel) screening in the high (low) density regime and 
the peak in the density dependence exactly occurs at 
the crossover between the degenerate and nondegenerate 
regimes. 

In contrast, electrons in most n-type zinc-blende ma- 
terials are easier to enter into the degenerate regime and 
the screening constant is usually small {n/kp < \/2) 
thanks to the small effective mass of the conduction band. 
Therefore, it is difficult to observe the nonmonotonic ef- 
fect of the screening during the electron spin relaxation. 
However, Jiang and Wu suggested that the screening 
from the holes in p-type materials can also give rise to ob- 
servable effects on the electron spin relaxation.— In that 
case, the EY mechanism is always irrelevant and the elec- 
tron spin relaxation is dominated by the DP mechanism 
in the low hole density regime. Since the spin relaxation 
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In the low density regime, the hole gas is in the nondegen- 
erate limit and the inhomogeneous broadening is small 
and insensitive to the density. Hence the contribution of 
the DP mechanism to the HSR is negligible. However, 
holes are in the degenerate limit for high density case and 
the inhomogeneous broadening increases rapidly [oc kp, 
see Eq. ^ also] with increasing hole density, which makes 
the DP mechanism markedly contribute to the HSR. By 
using Eq. (|15l) . one can easily obtain tdp/tey oc ?i^ 



for n/kp 3> v2 and tdp/tey 



oc n. 



for n/kp < \/2. It 



is obvious that tdp/tey decreases with increasing density 
in both cases. 
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FIG. 5: (Color online) (a) Density dependence of the total 
HSR time in p-type GaAs at 20 K. In the calculation, the 
RPA (red solid curve), Thomas-Fermi (green dotted curve) 
and Debye-Hiickel (blue dashed curve) screenings are applied, 
(b) HSR times due to the EY mechanism (curves with solid 
symbols) and the DP mechanism (curves with open symbols) 
as function of the hole density. 



time due to the DP mechanism is inversely proportional 
to the momentum scattering time in the strong scatter- 
ing limit, the hole density dependence of the electron spin 
relaxation time there is opposite to that of the HSR time 
we predict here. Specifically, they found that the elec- 
tron spin relaxation time due to the DP mechanism first 
increases and then decreases as the hole density increases 
in the nondegenerate regime. After the holes enter into 
the degenerate regime, the electron spin relaxation time 
increases again. 

From previous discussion, one notices that the HSR 
properties can be well interpreted by the EY mechanism 
which suggests that the EY mechanism is generally more 
important than the DP one. This can be clearly seen 
from Fig.[5lb), where the HSR times due to the EY and 
DP mechanisms are separately plotted. Interestingly, we 
find that the contribution of the DP mechanism can be 
comparable to that of the EY mechanism in high doping 
regime. To explain this behavior, we employ the relation 
in the strong scattering limit, tdp ^ l/[(|rip — D,1)Tp^] 
with (|r2p — n1) representing the ensemble average of the 
effective magnetic field (inhomogeneous broadening^^) . 
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FIG. 6: (Color online) (a) HSR times due to the EY (curves 
with solid symbols) and DP (curves with open symbols) mech- 
anisms as function of the temperature in p-type GaAs with 
Ui = Uh- The result without the hole-phonon scattering at 
Uh ~ 10^® cm~^ is plotted as dotted curve without symbol, 
(b) The ratio of the screening constant to the Fermi wave vec- 
tor as function of the temperature. The level k — ^plkp is 
illustrated by a solid line. 



2. Temperature dependence 

Similarly, one also finds the nonmonotonic tempera- 
ture dependence of the HSR time. In Fig.ElJa), we plot 
the results of four typical densities n^ — 10^^, 10^^, 10^^ 
and 4 x 10^^ cm~^. One notices that the temperature 



dependence of the HSR time is also qualitatively deter- 
mined by the EY mechanism (curves with solid sym- 
bols). The HSR time from this mechanism first de- 
creases and then increases with increasing temperature 
for Uh — 10^^ cm~^, and the minimum reaches around 
15 K. Since the hole-phonon scattering is rather weak in 
this regime, this feature just refiects the important role of 
the screening in the hole-impurity scattering. FigurelHUb) 
illustrates the ratio of the screening constant to the Fermi 
wave vector, where the crossover (k « y/2kp) occurs 
just around 15 K. Since holes are always nondegener- 
ate in this case (even at 5 K, see Fig.lll), the HSR time 
decreases (increases) with increasing temperature below 
(above) 15 K, according to Eq. p^ [Eq. pT|) ]. Moreover, 
we find that the HSR time decreases again when the tem- 
perature further increases, because of the enhancement of 
the hole-phonon scattering. To demonstrate this picture, 
we remove the hole-phonon scattering from the KSBEs 
and find that the HSR time due to the EY mechanism 
monotonically increases above 15 K, as expected (shown 
as dotted curve without symbol). For nt = 10^^ cm""^, 
only the peak is visible while the valley is absent as the 
curve with solid squares shown. The reason is that the 
screening is weak for this density [see Fig.[6Ub)]. How- 
ever, for Uh = 10^* and 4 x 10^^ cm^^, holes are in the 
degenerate regime at low temperature, hence the tem- 
perature dependence of the HSR time is determined by 
n/kp from Eq. (flS)) . At high temperature regime, the 
HSR time is also limited by the hole-phonon scattering 
for these densities. Another important information one 
can get from Fig.inja) is that the HSR time due to the 
DP mechanism is comparable to the one due to the EY 
mechanism around 80 K for Uh — 10^^ cm""^. For the 
high doping case with ti/i = 4 x lO""^^ cm~'^, the impor- 
tant role of the DP mechanism at low temperature can 
also be clearly seen [see also Fig.[5l[b)]. 
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FIG. 7: (Color online) Ratio of the HSR time due to the 
EY mechanism to that due to the DP mechanism, tey/tdp, 
as function of the doping density and temperature in p-type 
GaAs with Ui = Uh- The yellow dashed curves represent the 
borders with tey/tdp = 1/3 



ever, this effect is suppressed by the hole-phonon scat- 
tering at higher temperature, where the EY mechanism 
again is more efficient than the DP one. In the figure, 
we also plot the borders with tey/'^dp = 1/3 as yellow 
dashed curves. 

Finally, we should point out that the Bir-Aronov-Pikus 
mechanism^ is neglected in our computation even in the 
intrinsic case. The reason lies in the fact that the spin-flip 
process due to the electron-hole exchange interaction in 
intrinsic GaAs occurs in the time scale of nanoseconds^ 
therefore this mechanism is irrelevant to the ultrafast 
spin relaxation of the hole system. 



3. Relative contribution of the DP and EY mechanisms 



IV. CONCLUSION 



For thorough understanding of the relative contribu- 
tion of the DP and EY mechanisms, we plot the ratio 
of the HSR times due to the EY and DP mechanisms 
tey/tdp as function of the density and temperature in 
Fig-El It is clear to see that the HSR time is dominated 
by the EY mechanism at low temperature (below 40 K) 
upto 10^^ cm~^. This is because that holes mainly oc- 
cupy the states with small wave vectors and experience 
a small effective Dresselhaus field. In the high density 
regime, the increase of the inhomogeneous broadening 
makes the DP contribution significantly enhanced and 
even comparable to the EY one as discussed above. An- 
other interesting regime lies in the low density regime 
at moderate temperature, where the contribution of the 
DP mechanism can also be comparable to that of the EY 
mechanism. This originates from the reduction of the 
EY contribution due to the relatively weak hole-impurity 
scattering according to Eq. PT)) [see Fig.inja) also]. How- 



In conclusion, we have investigated the HSR from the 
fully microscopic KSBEs in intrinsic and p-type bulk 
GaAs. We analyze the valence-band structure by consid- 
ering the anisotropic property and the Dresselhaus spin- 
orbit coupling. We find that the degeneracy of the HH 
band and that of the LH band are both lifted by in- 
traband splittings. The DP mechanism associated with 
the intraband precessions and the EY mechanism due 
to the direct spin-flip scattering are then explicitly stud- 
ied, with the intraband coherence included. Our result 
of intrinsic GaAs shows good agreement with the exper- 
iment data at room temperature, where the EY mech- 
anism is demonstrated to be the dominant spin relax- 
ation mechanism due to the strong scattering process. 
We also show that the approach without the intraband 
coherence used in the literature is inadequate in account- 
ing for the HSR and the nonpolar hole-optical-phonon, 
missed in the previous theoretical work, is important to 
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the HSR. The temperature can markedly affect the HSR 
time by changing the strength of the hole-phonon scat- 
tering. At low temperature, we predict a valley due to 
the Coulomb scattering in the HSR time at the crossover 
between the degenerate and the nondegenerate regimes of 
the electrons. Inp-type GaAs, we find that the HSR time 
depends on density and temperature nonmonotonically. 
For the density dependence, the HSR time presents a 
peak at the crossover between the degenerate and nonde- 
generate regimes of the holes, resulting from the different 
features of the screening in the degenerate and nondegen- 
erate limits. In the nondegenerate regime, we predict a 
valley, from the competition between the screening con- 
stant and the momentum exchange, in the density depen- 
dence of the HSR time. We also find that the HSR time 
monotonically decreases with increasing density at high 
temperature, thanks to the small screening constant. For 
the temperature dependence, we predict a valley in the 
low temperature regime, which also reflects the role of 
the screening. In the high temperature regime, the hole- 
phonon scattering can markedly contribute to the HSR 
and make the HSR time decrease with increasing tem- 
perature. Moreover, we find that the contribution of the 
DP mechanism can be comparable with that of the EY 
one in the high density regime at low temperature and 
in the low density regime at moderate temperature, even 
though the EY mechanism is usually the major mecha- 
nism in the HSR. 



and those in Eq. ([5]) are given by Oj — a^ ^ h^ = 
-(61^')* and 4/' ==(ci/')*- 

The effective magnetic field in Eq. ^ can be written 



as 



^"^ = 



^l 



ni 



hJ{,e{\/ia2C2 - 62 + C2) - /i^Re(62C2) 

-I- ft.j^Im(\/3a2C2 -bl- cl), (A4) 

-/i2;lm(V3a2C2 - &2 + C2) + hzlm{b2C2) 



■ hyIie{V3a2C2 — b\ 



(A5) 



/ia;Re(v3a26i -I- 261C2) -I- /ij^Im(v3ai6i — 262C1) 
+ /i,(3a2-K|5ip-|cin/2, (A6) 



where the labels "/i/T' are neglected for short. 



Appendix B: Deformation potential matrix of 
optical phonons 



The 6x6 deformation potential matrix of optical 
phonons for valence bands reads^i 
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Appendix A: Detail of intraband splitting 

The coefficients in Eq. (|4]) can be expressed as 

i^E'^/i _ G)/[{E'''^ - Gf + \H\^ + |/p]^(Al) 



h/l 



,h/l _ 



h/l 



= H*/[{E''/' -Gf + lHl"^ 






(A2) 
(A3) 



with 



Tx^ = 



P: 



Aq 



a \ "Aq 02X2 



^Aq 









St 











\ 



-5i 



-^Sl -51 



¥x 



V25i 



¥1 



/ 



(Bl) 



2^A 



-V2<5f 



2 "A 



(B2) 



te3) 



where 5^ — 5\x ± i5\y and 5;^ = 5\z ■ do and a are optical 
deformation potential and lattice constant, respectively. 
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